Draft version May 21, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



oo 
O 
O 
(N 

>>■ 



(N 

6 



> 

OO 

o\ 
rn 

o 

oo 

o 



X 



SEARCH FOR BLUE COMPACT DWARF GALAXIES DURING QUIESCENCE 

J. Sanchez Almeida\ C. Munoz-Tunon^ R. Amori'n^, J. A. Aguerri^, 
R. Sanchez- Janssen\ and G. Tenorio-Tagle^ 

Draft version May 21, 2008 

ABSTRACT 

Blue Compact Dwarf (BCD) galaxies are metal poor systems going through a major starburst that 
cannot last for long. We have identified galaxies which may be BCDs during quiescence (QBCD), 
i.e., before the characteristic starburst sets in or when it has faded away. These QBCD galaxies 
are assumed to be like the BCD host galaxies. The SDSS/DR6 database provides '^21500 QBCD 
candidates. We also select from SDSS/DR6 a complete sample of BCD galaxies to serve as reference. 
The properties of these two galaxy sets have been computed and compared. The QBCD candidates 
are thirty times more abundant than the BCDs, with their luminosity functions being very similar 
except for the scaling factor, and the expected luminosity dimming associated with the end of the 
starburst. QBCDs are redder than BCDs, and they have larger H ii region based oxygen abundance. 
QBCDs also have lower surface brightness. The BCD candidates turn out to be the QBCD candidates 
with the largest specific star formation rate (actually, with the largest Ha equivalent width). One 
out of each three dwarf galaxies in the local universe may be a QBCD. The properties of the selected 
BCDs and QBCDs are consistent with a single sequence in galactic evolution, with the quiescent phase 
lasting thirty times longer than the starburst phase. The resulting time-averaged star formation rate 
is low enough to allow this cadence of BCD - QBCD phases during the Hubble time. 
Subject headings: Galaxies: starburst - Galaxies: evolution - Galaxies: dwarf - Galaxies: luminosity 
function - Galaxies: abundances 



1. INTRODUCTION 

Blue Compact Dwarf (BCD) galaxies'^ are met al poor 
systems undergoing vigorous star formation fe.g.. lThuanl 
119911: iGil de Paz et al.ll2003[) . W ith record-bre aking low 
metallicities among galaxies fe.g.. lKunth fc Os thn 2000), 
their observed colors and spectra indicate essentially 
newborn starbursts, with r nean ages of a few My (e.g., 
iMas-Hesse fc KunthI Il99l iThuanI |1991D . This combi- 
nation of factors (chemically unevolved systems, with 
oversized starbursts that cannot last for long) lead to 
conjecturing that they are pristine ga laxies undergoing 
star formation for the very first time (jSargent fc Searlel 
119701 ). This original view has now been reformulated 
so that BCD galaxies are chemically primitive objects 
which we come across during short intense bursts. The 
starbu rst phases are inte rleaved by long periods of quies- 
cence (jSearle et al.lll973l) . The change is based on several 
observational evidences. Most BCDs galaxies are known 
to have a red lo w surface brightness component (e.g., 
ICaon et al.l I2005L and references therein), which should 
exist before the starburst, and which should survive the 
BCD phase. We are also allowed to resolve some indi- 
vidual stars in the nearest BCD galaxies, and the pres- 
ence of RGB stars indicates an underlying stellar popu- 
lation much older than the lifetime of a starburst (e.g., 
lAloisi et al.l[2007l : iCorbin et al.l [200?). 

Even if they are not pristine, BCDs are galaxies with 
the lowest metallicities and may therefore be showing 



first stages in star formation from primordial gas. It is 
not yet clear which objects eventually change to glow 
as BCD galaxies (i.e., which galaxies are BCDs dur- 
ing quiescence). Identifying them seems to be a nec- 
essary first step to answer the question of why a galaxy 
experiences a BCD phase. What is the nature of the 
host galaxy in which a BCD starburst takes place? 
Is it a normal dwarf galaxy? What is triggering the 
starburst? What is left after the starburst? Evolu- 
tionary connections between dwarf elliptical galaxies, 
dwarf irregular galaxies, low surface brightness galax- 
ies, and BCDs have b een both proposed and criticized 
in th e literature (e.g., Searle & Sargent 1972'; Silk et aL| 
19871 [bavies fc PhiUipps 198 , 8; P apadcros et al. 
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Here we take a new approach to investigate the linage 
of the BCDs. Rather than directly comparing the prop- 
erties of the BCDs (or their host galaxies) with proper- 
ties of known galaxy types, we attempt a blind search 
for galaxies that look like precursors or left-overs of the 
BCD phase, paying minor attention to the galaxy class. 
We have tried to find field galaxies with the properties of 
the BCD host, but without the BCD starburst. Isolated 
galaxies are preferred to minimize the role of mergers 
and harassment on the galactic properties and evolution. 
These galaxies will be denoted in the paper as Quies- 
cent BCD galaxies or QBCD. If there are enough such 
galaxies, the BCD may be just a particularly conspicu- 
ous phase in the life of otherwise mean dwarf galaxies. If 

.WyjrbBt@ittHeS^pt®fihki*iEftix]^^ should conclude that the 
ULJJJ phase enas up with the host, which is physically 
unlikely. Aiming at shedding light on the BCD evolu- 
tion, we have tried to answer two specific questions: (1) 
are there galaxies like the BCD hosts without the con- 
spicuous starburst observed in BCDs? (2) If the answer 
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were yes (as it turns out to be), what are the physical 
properties of these (putative) BCD galaxies during qui- 
escence? 

The different sections of this paper describe individ- 
ual steps carried out to answer the two questions posed 
above. The search for QBCD galaxies like the BCD 
host galaxies is described in § 12.11 Luminosities, sur- 
face brightnesses, and c olors character i stic o f BCD host 
galaxies are taken from I Amorin et al.l (|2007l [2008). We 
search the SDSS/DR6 database for candidates, which is 
ideal for the kind of exploratory statistical study we aim 
at. The sample of QBCD galaxies has to be interpreted 
in terms of a sample of BCD galaxies. In order to min- 
imize systematic errors, the control sample of BCDs is 
also derived from SDSS/DR6 using the same techniques 
employed for QBCD selection (§ 12. 2p . We compare the 
Luminosity Functions (LFs) of the two samples in § [31 
which requires computing how the BCD galaxy luminos- 
ity changes to become a QBCD galaxy when the star- 
burst fades away (333]) • Since the SDSS catalog is mag- 
nitude limited, the properties that we derive are biased 
towards the brightest galaxies (Malmquist bias). This 
effect is corrected to derive the intrinsic properties of 
the two sets (§ S]). The implications of our findings are 
discussed in § [5] and § [6l We have developed our own 
software to compute LFs. For the sake of clarity and 
future reference, details are given in App. |^ A Hubble 
constant Ho — 70 km s~^ Mpc~^ is used throughout the 
paper. 

2. DATA SETS 

The search has been carried out using the Sloan 
Digital Sky Survey (SDSS) dataset, which is 
both convenient and powerful. We use the latest 
data release, DR6, whose spectroscopic view cov- 
ers 7425 deg^ and contains ~ 7.9 x 10^ g alaxie s 
(jAdelman-McCarthv k the SDSS collaboration! I2007D . 
The main characterist ics of the SDSS are de scribed in 
the extensive paper bv lStoughton et al.l (|2002f) . but they 
are also gathered in the comprehensive searchable SDSS 
website^. We access the database through the so-called 
CasJob entry, which allows direct and flexible SQL 
searches, such as those needed to select only isolated 
galaxies. Absolute magnitudes are computed from 
relative magnitudes and redshifts. All magnitude s have 
been corrected for galactic extinction (Sch lcgel et al.l 
I1998D . No ivT— correction has been applied since we 
are dealing with nearby galaxies with redshifts < 0.05, 
and so, with Hubble flow induced bandshifts much 
smaller tha n the relative bandwidth of the color filters 
(jFukugita et al. 1996). In a consistency test using SDSS 
spectra representative of QBCDs and BCDs, we have 
estimated that the iiT— corrections are smaller than 0.01 
magnitudes (QBCDs) and 0.1 magnitudes (BCDs). 

2.1. Selection of galaxies like the hosts of BCDs 

The properties of the QBCD galaxies have been 
taken from the sample of 28 BCD galaxies in 
ICairos et al.l (|2001D. whose hosts h ave been character- 
ized bv lAmorfn et all (|2007l |2008[ ). They represent a 
large set of host galaxy 
from 2-dimensional fits. 
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Fig. 1. — Scatter plots with magnitudes, colors and concen- 
tr ation indexes ch aracterizing the hosts of BCD galaxies according 
to lAmorm et all EOOTI . 1200a) . (a) Color vs absolute magnitude, 
(b) Redest SDSS magnitude z (9097 A) vs visible magnitude g 
(4825 A), (c) Sersic index (= concentration) vs magnitude, (d) 
Sersic index vs color. The plus symbols inserted into the boxes 
mark the gold subsample, i.e., those galaxies which arc believed to 
be characterized best. 



Sersic profiles^ to the outskirts of each BCD galaxy once 
the bright central blue component has been masked out. 
From fits in various bandpasses, they derive luminosi- 
ties, colors and Sersic indexes. We use them to guide the 
SDSS search. (The impact on the selection of using this 
particular characterization rather than other alternatives 
is discussed in § 15.41 ) The actual properties of these BCD 
host galaxies are shown in Fig. [Tj The SDSS photometric 
col ors uqriz are not the standard UBVRI system used 
by lAmorm et al.l Therefore, in order to produce Fig. [1] 
the original de- reddened colors have been translated to 
SDSS magnitudes using the recipe in iSmith et al.l (|200l 
Table 7). This photometric transformation is suited for 
stars but, according to the SDSS website, the equations 
also provide reliable results for galaxies without strong 
emission lines, as we expect the BCD hosts to be. 

Note that the absolute magnitude, the color, and the 
Sersic index of each BCD host do not seem to be corre- 
lated with each other (cf. Figs. [T^, [TJ:, and[TJi). This 
fact allows us to assume that the magnitude, color and 
concentration of each host galaxy are independent, and 
we will look for galaxies whose properties span the full 
range of observed values. 

Before we can proceed, there is an important obser- 
vational property of the BCD hosts tha t needs to be 
ment ioned. (For a detailled discussion, see lAmorfn et al.l 
l2008l ) BCD galaxies and their host galaxies tend to 
have similar luminosities, so tha t the brighter the BCD 
galaxy the brighter t he host (see iPapaderos et"al]ll996al : 
iMarlowe et alj [i999[ ). Actually, there seems to be a 
simple relationship between the magnitude of the host, 
Mhost, and the magnitude of the BCD galaxy, AfecD, 



Afhost ^ 0.5 + Mbcd. 



(1) 



' |http : //mv . sdss . org/dr6 1 



^ The magnitude m varies with the radial distance from the 
center of the galaxy p as m{p) = m(0) — bn[p/ Re]^^" , with b„ 
chosen so that half of the galaxy light is enclosed within p < Re 
(see, e.g.. [Ciottiiil99lA . The so-called Sersic index n controls the 
shape of the profile. 
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Fig. 2. — Absolute magnitude of the low surface brightness host 
of the BCD vs magnitude of the BCD, including the host. Each 
symbol corresponds to a Johnson's bandpass as labeled in the inset. 
There is a tight correlation between the magnitude of the BCD and 
the magnitude of the host galaxy. A linear least squares fit based 
on B magnitudes yields a slope of 0.991 it 0.054, and an offset of 
0.61 ± 0.11 (the dashed line). Within error bars, the results are 
identical when using the full data set for fitting (the dotted line), 
indicating a relationship i ndepe ndent of the color. The gold set 
defined by lAmon'n et al.l 112008) also yields the same result (not 
shown) . 



TABLE 1 

Criteria used to select BCD host-like galaxies (i.e., 
QBCDs) FROM THE SDSS/RD6 database. 



# 


Criterion 


Implementation 


(1) 


Colors 


-0.5 < Mg- Mr < 1.1 mag 


(2) 


Concentration indexes 


1.8 < R90/R50 < 3.1 


(3) 


Magnitudes 


-18.5 < Mg < -13.5 mag 


(4) 


Be isolated 


no bright galaxy within 10 


(5) 


Get rid of proper motion 


redshift > 0.003 


induced redshifts 





Bright means brighter than 3 magnitudes fainter than the selected 
galaxy 

Figure [5] presents the scatter plot bet ween the ob- 
served magnitude of the BCDs analyzed bv lAmorm et al.l 
()2008f ). and the magnitude of their host galaxy. There 
is a linear relationship between the two observed mag- 
nitudes as parameterized by equation ([T]) ~ the slope is 
very close to one, and the scatter remains rather small 
(less than 0.5 mag). The relationship is independent of 
the color and, therefore, it will remain like equation ([1]) 
for the SDSS photometric system. 

Table [T] lists the actual constraints employed in our 
SDSS/DR6 search for QBCD galaxies like the BCD host 
galaxies. The range of colors, 



0.5 < A/„ - Mr < 1.1, 



(2) 



has been taken directly from Fig. [T^. The range of ab- 
solute magnitudes, 

- 18.5 < Mg < -13.5, (3) 

requires a more elaborated explanation. The upper 
bound corresponds to the magnitude of the faintest BCD 
host (Fig. [1^). The lower limit, however, is inherited 
from the lower limit of the BCD galaxies selected in § 12.21 



This is a control set used for reference and, according to 
arguments to be given in § 12.21 SDSS BCD galaxy can- 



didates are chosen to be fainter than Mg ~ —19. There- 
fore, for consistency with equation ([T]), we impose the 
lower limit given in equation ([3|). There is an additional 
detail conce rning the magnitudes in use. I Amor m et al.l 
(|2007l l2008f ) integrate the Sersic profile to infinity to es- 
timate magnitudes, whereas the SDSS catalog provides 
Petrosian magnitudes, where the galactic light is inte- 
grated onl y up to a certain dista nce from the g alactic 
center (see lstoughton et al.|[200l ). iGraham et alI(|2005D 
show how the difference between the Petrosian magni- 
tude and the total magnitude differs by less than 0.2 mag- 
nitudes for Sersic indexes n < 4. The hosts have Sersic 
indexes smaller than this limit (Fig. [T};) and, therefore, 
we neglect the difference. 

The brightness profile of the B CD host is parameter- 
ized in lAmorm et al.l (|2007l I2OO8I ) by means of the Sersic 
index n. According to Fig. [H;, 

0.5 < 7i< 2.9. (4) 

This constraint is translated into a limit on the so-called 
concentration index Rqq/R^q, which can be readily ob- 
tained from the SDSS catalog since it provides the radii 
containing 50% of Petrosian flux, Rr,^, and 90% of the 
Petrosian flux, Rqq. According to iGraham et al.l (|2005l 
Table 1), the limits given in equation (j4]) imply, 



1.8 < i?go/ R50 < 3.1, 



(5) 



which is the constraint used in our search (Table [U 2nd 
row). This constraint has been applied both to g and to 
r magnitudes. 

As explained in the introduction, we seek isolated 
galaxies to minimize the effect of mergers and harass- 
ment on the galaxy properties. The criterio n for finding 
isolated galaxies is inspired by the work of lAllam et al.l 
()2005f ). We stipulate that the selected galaxies have no 
bright companion within 10 R^o- (Note that R50 is ap- 
proximately the effective radius, i.e., the radius contain- 
ing half of the galaxy luminous flux.) Companions are 
neglected if they are at least three magnitudes fainter 
than the candidate. These two constraints are set in the 
color filter g (see Table [1] 4th row). 

Since absolute magnitudes are computed from red- 
shifts, we ask the redshifts to be large enough to minimize 
the proper- motion induced Doppler shifts. The thresh- 
old redshift corresponds to a distance of 13 Mpc and a 
velocity of 900 km s"^ (see TablelU 5th row). 

Applying the criteria in Table[T]to the SDSS/DR6, one 
retrieves 21493 galaxies. Their mean redshift is 0.030, 
with a standard deviation of 0.014. The galaxies are illus- 
trated in Fig.[3l top row, where we include color images of 
four randomly chosen QBCD candidates. Figure 2] shows 
various histograms corresponding to the observed proper- 
ties of this set of galaxies (the solid lines). Means, stan- 
dard deviations and modes corresponding to these his- 
tograms are listed in Table [2l The determination of the 
concentration indexes of small galaxies may be problem- 
atical. They are based on determining ga laxy sizes and, 
therefore, they are affected by seeing (see iBlanton et al.l 
l2003af ). We select a subset among the QBCD candidates 
to test whether the main properties of the candidates de- 
pend on the spatial resolution. To be on the safe side, 
we consider candidates having i?5o > 4", i.e., larger than 
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Fig. 3. — Set of randomly chosen images of QBCD candidates (top) and BCD candidates (bottom). We are using the color loolc-up 
table of SDSS. Each galaxy includes a horizontal scale corresponding to 5". The crosses centered in the galaxies point out north, south, 
east and west directions. 



the typical SDSS seeing (< 1.5"). The corresponding 
histograms are also shown in Fig. |4] (the dashed lines), 
with the means, standard deviations and modes included 
in Table O Except for the surface brightness, there are 
no systematic differences between the full set and the 
subset of large QBCD galaxies. Moreover, the decrease 
of surface brightness for large galaxies has nothing to 
do with poor seeing. It is a bias imposed by our lower 
limit in absolute magnitude (Table [U 3rd row). Unless 
the large galaxies are also low surface brightness, they 
are much too bright to satisfy our selection criterion. In 
short, poor seeing does not seem to bias our selection in 
any obvious way. 

2.2. Selection of BCD galaxies 



The galaxies in §[27T]are selected to be quiescent BCDs. 
It is clear that this conjecture and the galaxy properties 
must be examined in terms of the properties of the BCD 
galaxies. This is particularly true to answer the first ba- 
sic question of whether the number density of QBCD can- 
didates suffices to account for the existing BCDs. Num- 
ber densities are best characterized as LFs but, to the 
best of our knowledge, there is no BCD LF in the litera- 
ture. Moreover, even if such a LF existed, it would have 
been produced with a number of (unknown) biases dif- 
ferent from those involved in our SDSS selection. There- 
fore, we found it necessary to extract from SDSS a set 
of BCD galaxies, that is to say, to construct a reference 
BCD sample with the biases and problems of the QBCD 
galaxies we want to compare to. This section describes 
such selection of BCD candidates. 

The general criteria have been taken form iMalmbergl 

(l2005f). but th ey are implemented according to 

IGil de Paz et al.l (j2003f ). A BCD galaxy should have 
the following properties, (1) be blue enough, which con- 
strains the colors, (2) be compact, which limits the sur- 
face brightness, (3) be dwarf which sets a lower limit 



to the magnitude, (4) have a large star formation rate 
(SFR), which implies having enough H ii regions and, 
therefore, enough Ha emission, (5) be metal poor, and 
(6) not to be confused with AGNs. These general criteria 
have been specified as detailled in Table [3l The numer- 
ical values correspon ding to constraints (1), (2) and (3) 
have been taken from lGil de Paz et al.l ()2003l ). who point 
out three photometric criteria for a galaxy to be a BCD, 
{i^b) — (ij-r) < 1 niag arcsec"^, 

{fis) < 22 mag arcsec"^, (6) 
Mk > -21mag. 

Using the transformatio n between Johnson 's and SDSS 
photometri c systems in ISmith et aH ()2002f ). plus equa- 
tion (1) in IGil de Paz et al.l l20oW . these three condi- 
tions become the constraints (1), (2) and (3) in Tabled 
The mean surface brightness (fi) has been computed as 
the magnitude of the average luminosity within i?5o, i.e., 
(/i) = 2.51og(27r) + m + 51ogi?5o, with m the appar- 
ent magnitude (see, e.g.. lBlanton et aI1l2001l § 2.3). The 
transformation between photometric systems is suited for 
stars, rather than for galaxies with emission lines like the 
BCDs. However, the approximation suffices because the 
actual limits in equation ([6]) are estimative, and the effect 
of including lines modifies the BCD magnitudes by 0.1 
mag or less. This effect has been estimated using SDSS 
spectra representative of BCD galaxies. 

The actual implementation of criteria (4) and (5) are 
mere educated guesses that try not to be too restrictive. 
For example, we consider metallicities smaller than 1/3 
the so lar value (Table |3l item 5), when iKunth fc OstlinI 
(|2000( ) bound the BCDs metallicities between 1/10 and 
1/50. The (oxygen) metallicity of the galaxies (O/H) has 
been estimated using the so-called N2 method, based on 
the equivalent widths of the emission lines [A^//]A6583 
and H„, 



12 + log(0/H) = 9.12 + 0.73 N2, 



(7) 
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Fig. 4. — Histograms of the observed magnitudes and colors corresponding to our selection of QBCD galaxies (the solid lines), and BCD 
candidates (the dotted lines). We also consider a sub-set of QBCD galaxies with large apparent size (Rso > 4"). Their properties are 
shown using dashed lines. The labels in the ordinate axes specify the parameter that is represented - from top to bottom and left to right, 
absolute z magnitude Afz, absolute g magnitude Mg, surface brightness color {fig) — {p.r), surface brightness (fJ.g), oxygen-to-hydrogen 
abundance in the usual scale 12 -|- log(0/H), and Ha equivalent width. The histogram of metallicities also includes the solar value for 
reference (12 4- log[0/H]Q ~ 8.91: IKunth fc OstliD|[2000t) . 



TABLE 2 

Statistical properties of various galaxy sets. 



Galaxy 






Mg 




(fig) - (fir) 






12+log(0/H) 


Ha 




Set 










[mag arcsec 




[mag arcsec 








[A] 






av'^itstd'^ 


md<i 


avistd 


md 


avitstd 


md 


avitstd 


md 


avitstd 


md 


avitstd 


md 


Observed QBCD 


-18.5±1.1 


-18.8 


-17.6±0.8 


-18.3 


0.50±0.24 


0.46 


21.8±0.88 


22.0 


8.61it0.30 


8.60 


28.1it49.2 


1.9 


Large QBCD 


-18.0±1.3 


-18.6 


-17.3±1.0 


-18.1 


0.47±0.16 


0.46 


22.7±0.67 


22.9 


8.61it0.26 


8.60 


19.9it35.5 


9.9 


Observed BCD 


-17.9±1.1 


-18.3 


-17.5±1.1 


-19.0 


0.04±0.23 


0.15 


20.8±0.54 


21.0 


8.24it0.16 


8.35 


172. ±174. 


58. 


Restored QBCD 


-16.2±1.6 


-14.9 


-15.5±1.4 


-14.2 


0.37±0.23 


0.38 


22.2±0.95 


22.6 


8.43it0.31 


8.40 


33.1it50.4 


1.9 


Restored BCD 


-16.2±1.5 


-16.3 


-15.8±1.4 


-15.1 


0.04±0.22 


0.15 


21.0±0.59 


21.7 


8.12it0.20 


8.14 


121. di 64. 


58. 



^ Ha equivalent width Average'^ Standard deviation Mode 



with 



N2 log 



[iV//]A6583 

Ha 



see iShi et all ()2005D and references therein. Note that 
this metallicity characterizes the properties of the H ii 
regions excited by the starburst, which represents only a 
very small fraction of the galactic gas (see § O . Keeping 
the appropriate caveats in mind, the N2 method suffices 
for the elementary estimate we are interested in. 



Active Galactic Nuclei (AGNs) can be misidentified as 
star-forming galaxies since they are high surface bright- 
ness blue galaxies with emission lines. AGNs introduce 
spurious BCD candidates. In practice, however, this con- 
tamination is negligible, and we do not decontaminate 
for the presence of AGNs (item 6 in Table [3]). Our ap- 
pr oach can be justified usin g the criterion to be an AGN 
bv iKauffmann et al.l ()2003l ). An emission line galaxy is 
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TABLE 3 

Criteria used to select BCD galaxies from the SDSS/DR6 databas]ec 



# 



Criterion 



Implementation 



(1) Be blue enough 

(2) Be compact 

(3) Be dwarf 

(4) Having large SFR 

(5) Be metal poor 

(6) Not to be confused with AGNs 

(7) Be isolated 

(8) Get rid of proper motion 
induced redshift 



(/jg) — (fir) < 0.43 mag arcsec ^ 
{fig) < 21.83 ~0.47{(Mg> - ifir)) ngi! 
Mg > -19.12 + 1.72{Mg - Mr) m^ 
Ha Equivalent Width > 50 A 
12 + log{0/H) < 8.43 {= 1/3 0) g 
neglect AGN contamination 
no bright galaxy within 10 -Rso^ 
redshifts > 0.003 



nagnitudes faint c 



L the selected galaxy. 



an AGN if 

, [O///]A5007 ^ ^ 
log^ -r; > 1-3- 



13 



0.61/ log 



[iV//]A6583 



0.05 



(8) 

We sel ected all SPSS BCD c andidates chosen accord- 
ing to iGil de Paz et al.1 (|2003D criteria (items 1 , 2 and 
3 in Table [3]), and having the spectral line informa- 
tion required to apply the test (the hnes [O///]A5007, 
[A''J/]A6583, Hon and Hf^). Figure [5] shows the scatter 
plot of the two indexes involved in equation ([8]) . Points 
above the solid line are AGNs according to the criterion 
in equation ([S]). There are very few AGNs and, what is 
even more important, most of them stay to the right of 
the vertical dashed line, i.e., 

, [iV//]A6583 
N2 = log ^ ^ > -0.95, (9) 

which corresponds to the constraint on the N2 index im- 
posed when the 0/H abundance estimate is based on 
the N2 index, and the O /H is constrained as we do (item 
5 in Table [3]). Our BCD candidates are to the left of 
this line and therefore the low 0/H abundance constraint 
automatically removes most of the AGN contamination, 
which explains our approach. In order to quantify the 
residual contamination, let us mention that only 0.4% of 
the low 0/H abundance BCD candidates in Fig. [5] are 
also AGN candidates. 

We select isolated BCD candidates (criterion 7 in Ta- 
ble [3]) because of consistency with the criteria used to 
search for BCD host galaxy candidates. As in the case 
of QBCDs, we stipulate that the galaxies have no bright 
companion within 10 i?5o, where i?5o is the radius in- 
cluding 50% of the Petrosian flux. The companions are 
not bright enough to perturb the galaxy if they are three 
magnitudes fainter than the galaxy. 

Applying the criteria in Table[3]to the SDSS/DR6, one 
gets 1609 BCD galaxies. The dotted lines in Fig.[l]repre- 
sent the distribution of physical properties of this set of 
candidates. Means, standard deviations, and modes are 
also included in Table [H Four randomly chosen BCD 
candidates are shown in Fig. [3l bottom. The main red- 
shift of the BCD candidates is 0.032, with a standard 
deviation of 0.019. 

After carrying out the selection described above, we 
figured out that most of the BCD sample is also part of 
the QBCD sample: 1198 galaxies are shared by the two 
sets. They represent 74.5% of the BCDs and 5.6% of 
the QBCDs. The typical properties of the sets are very 
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Fig. 5. — Diagnostic plot used to identify AGNs contaminating 
the BCD sample. The points above the thick solid line correspond 
to AGNs candidates. The vertical dashed line separates the low O 
abundance region (to the left) and the high O region (to the right). 
We only use low O candidates, where the AGN contamination is 
minimum. The two indexes are defined in the main text. 

different (Table [2 and Fig. 2]) , however there is overlap- 
ping between the two populations. Roughly speaking, 
the BCD sample represents the fraction QBCD galaxies 
having the largest H^ equivalent width (see Fig. IH bot- 
tom right, where the dotted line and the solid line agree 
ioxHa EW> 100 A). Since the Hg lum inosity is a proxy 
for star formation (e.g. ' KennicuttllT998f ) . the BCD sample 
seems to be the QBCDs having the largest specific Star 
Formation Rate (SFR), i.e., the largest SFR per unit of 
luminosity. The overlapping is consistent with the two 
galaxy sets being part of a single continuous sequence, 
the most active QBCD galaxies still being identified as 
BCD galaxies. 

3. LUMINOSITY FUNCTIONS 

The LF, <I>(Af), is defined as the number of galax- 
ies with absolute magnitude M per unit volume and 
unit magnitude. We will use it to quantify the num- 
ber of QBCD candidates per BCD galaxy existing in the 
nearby universe. LFs are computed using a maximum 
likelyhood procedure s imilar to the methods described 
in the literature (e.g.. Efstathiou et al."l988': iLin et al.l 
[1996; Takcuchi ct al. 2000; Blanton et al. 200 11). We de- 
veloped and tested our own code for training purposes, 
and its characteristics and similitude with existing meth- 
ods are described in App. El Cubic splines are used to 
parameterize the LF shape, and the bes t fit is retrieved 
maximizing the likelyhood in equation (jA9|), a task re- 
alized with the usual Powell algorithm (e.g. iPress et al.l 
Il988f) . Errors bars are assigned by bootstrapping. 

When the procedure is applied to the QBCD galaxies 
described in § 12.11 one finds the LF represented in Fig.[6K 
(the solid line) . This LF is similar to that obtained with 
the 1/Umax method working on the same data sets (see 
App.[A|. The formal errors deduced from bootstrapping 
are very small (the shaded area around the solid line in 
Fig-Hfc-)- We use 50 bootstrap re-samples, but the error 
estimate is not very sensitive to this parameter. The LF 
is mostly sensitive to the apparent magnitude limit of 
the QBCD galaxy set. We are using a single apparent 
magnitude limit for the full dataset, and changing this 



BCD galaxies during quiescence 



7 



2.5e-02 QBCD galax Mpc-^ Observed 



10000 




Fig. 6. — (a) g color LFs of the QBCD galaxy sample (the solid line), and the BCD galaxy sample (the dashed line). The shaded area 
shows our statistical error estimate. The dotted-dashed line represents the LF to be expected when the starburst of the BCD galaxi e s fade s 
away revealing the underlying host galaxies. The dotted line corresponds to the LF of low luminosity galaxies bv iBlanton et ahl II2005I ). 
and it is shown here for reference, (b) Observed apparent magnitudes g, with the vertical lines corresponding to the magnitude limits 
used to derive the LFs. As the inset in (a) indicates, the solid line and the dashed line represent the QBCD sample and the BCD sample, 
respectively, (c) Scatter plot of the absolute magnitude of the QBCD galaxies vs the redshift. The solid line shows the boundary to be 
expected if the sample were limited in magnitude with g < 18.3. 



limit modifies the overall normalization - given a number 
of observed galaxies, the deeper the magnitude limit the 
lower the inferred number density of galaxies. We take 
for the limit, 

g < 18.3. (10) 

This selection is consistent with the magnitude 
limit of the main SPSS galaxy sample (r < 17.8 , 
lAdelman-McCarthv fc the SPSS collaboration! I2007D . 
keeping in mind that the QBCP galaxies are somewhat 
red with g — r ~ 0.5 (see Table [5] and Fig. 2]). However, 
we obtain the limit in equation (jlOp from the scatter plot 
of the QBCP galaxy absolute magnitude vs the redshift 
shown in Fig.[6j:. By trial and error, we modify the curve 
representing the boundary to be expected if the whole 
sample were limited with a single apparent magnitude. 
The best match is shown as the solid line in Fig. [6}: and, 
except for the range of large luminosities, the data set fits 
in well the expected behavior. Galaxies whose apparent 
magnitudes exceed the threshold in equation PU)) are not 
used to compute the LF - a histogram of the observed 
apparent magnitudes is shown in Fig.O) (the solid line). 
The normalization of the LF depends on the fraction of 
sky covered by survey, which we take to be the coverage 



of the SPSS/DR6 spectroscopic catalog (7425 d eg^; 
lAdelman-McCarthv fc the SPSS collaborationll2007[) . 

Figure [6^ also shows the LF corresponding to the BCD 
galaxies selected in 12.21 (the dashed line). The number 
density of BCP candidates {uq bcd — 9.2 10^** Mpc~'^) is 
much smaller than the number density of QBCP galaxies 
{no QBCD ^ 2.5 10^2 Mpc-3), giving, 

noQBCD^27. (11) 

'^O BCD 

In this case the limit magnitude of the sample has been 
set to g < 18.5, which we determine with the same pro- 
cedure used for the QBCD galaxies. The need for a 
different apparent magnitude limit can be inferred from 
the histograms in Fig. [6)3. Note how the histogram for 
BCD galaxies (the dashed line) does not drop off for large 
magnitudes as abruptly as the histogram for the QBCD 
galaxies (the solid line). 

The calculus of LFs has been repeated using z magni- 
tudes, i.e., the redder SDSS color filter at 9097 A. This 
bandpass is less sensitive to (blue) starbursts, enhancing 
the contribution of old stellar populations. The results 
are represented in Fig. [71 They are similar to those ob- 
tained from g, except that the LFs are shifted by one 



8 



Sanchez Almeida et al. 




Fig. 7. — This figure is identical to Fig. [6l except that the LFs, the histograms, and the scatter plot refer to z, i.e., the redder color 
among the SDSS bandpasses. See Fig. |6] caption for details. 



magnitude. The ratio between the number of QBCD 
galaxies and BCD galaxies remains as in equation (fTT|) - 
24 in this particular case. 

Figures [6^ and [7^ include the LFs for extremely low 
luminosity galaxies worked out by iBlanton et al.l ()2005l ) 
(the dotted hues). We used them as a reference, since 
they include all low redshift galaxies in the SDSS spec- 
troscopic catalog and, therefore, they provide LFs for the 
local universe affected by the same kind of bias as our 
galaxy selection. Using these LFs as reference, we find 
the QBCD galaxies to be rather numerous. One out of 
each three dwarf galaxies is a QBCD candidates (cf. the 
solid lines and the dotted lines in Figs. [5^ and [7^). Sim- 
ilarly, one out of each ninety dwarf galaxies is a BCD 
candidate^ (cf. the dashed lines and the dotted lines in 
FigsEJa and Ell). 

Our LFs are not corrected for the incompleteness of 
the SDSS catalog at low surface brightne ss. According 
to the detained simulation carried out bv IBlanton et al.l 
(|2005l Fig. 3), the incompleteness is not significant for 
galaxies with (/ig) < 23.5. Most of our QBCD candi- 
dates are within this bound (Fig.[4l bottom left), and it 
does not affect the BCD selection at all (Table The 
QBCD LF would be affected at its faint end, but even 
at Mg ~ — 14 the expected correction is not larger than 

® In a reply to J. Young. [ThuanI 1119911 ) estimates the same ratio 
using very different arguments. 



a factor of two (jBlanton et al.l [20051 . Fig. 6). This un- 
certainty does not modify the conclusions of the paper, 
which are mostly qualitative. However, the number of 
QBCD galaxies worked out in the paper may be under- 
estimated by a factor of order one. 

3.1. From LF o/BCDs to LF of BCD hosts 

The question arises as to whether the LFs for QBCD 
galaxies and BCD candidates are consistent. They seem 
to be, according to the following heuristic procedure to 
estimate the LF for the hosts of the BCDs. Given a LF 
of BCDs, $BCD, one can work out the LF to be expected 
once the starburst dies out, $qbcd- By switching the 
star formation on and off, BCD galaxies turn into QBCD 
galaxies, and vice versa. These two evolutionary phases 
have different timescales; tbcd for the BCD phase, and 
tqbcd for the QBCD phase. In a stationary state the 
number of galaxies turning from BCD to QBCD must be 
balanced by the galaxies going from QBCD to BCD, 

$qbcd(M)^^ ^ $bcd(M')^^, (12) 
T'QBCD Tbcd 

where BCD galaxies with magnitude in the interval 
Af ± AAf /2 become QBCD galaxies with magnitudes 
M± AM/2. Combining equation with the empirical 
relation between the magnitudes of the host galaxy and 
the BCD (equation [1]), one has the recipe to infer the 
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LF of the host galaxies to be expected from the LF of 
the BCD galaxies, 



$qbcd(M) 



TQBCD 
TBCD 



*bcd(M-0.5). 



(13) 



The scaling factor tqbcd/'''bcd is just the ratio between 
the number densities of BCDs, hq bcd, and QBCDs, 

no QBCD, 

'''QBCD no QBCD 



tbcd 



no BCD 



(14) 



since $qbcd and $bcd are normalized to no qbcd and 
no BCD, respectively. 

Using equations and (fT4)) . we evolve the LF of 
BCD galaxies in Figs. [S^ and [7^ to obtain the LFs of 
their host galaxies. They are included in the same fig- 
ures as the dotted-dashed lines. The scaling has 
been taken as the ratio between the number densities 
of QBCD galaxies and BCD galaxies (i.e., the ratio in 
equation |11] for the g filter, and the corresponding fig- 
ure for the z filter). The similarities between the evolved 
LF (the dotted-dashed lines) and the LF computed from 
the QBCD data sets (the solid lines) are quite strik- 
ing. We interpret this agreement as an indication of self- 
consistency between the two sets of galaxies chosen in 
§ [21 BCD galaxies and QBCD galaxies can be different 
phases in the life of a dwarf galaxy, with the lifetime in 
the BCD phase some thirty times shorter than the life- 
time in the QBCD phase (equations [11] and [14]). 

The BCD starbursts are very young since they still 
conserve massive stars (see §[T]). If tbcd — 10 My then 
equations ([T4]l and ([TT]) imply tqbcd — 0.27 Gy. Con- 
sequently, each QBCD galaxy may undergo as many as 
30-40 star formation episodes during the timespan where 
stars can be formed in dwarf galaxies (>10 Gy; e.g., 
iKunth fc Osnh3[200(Th . This issue is discused in § [5] 

4. PROPERTIES OF QBCD GALAXIES AND BCD 
GALAXIES 

The histograms of galaxy properties discussed in § [S] 
and shown in Fig. [H describe the properties of the ob- 
served galaxies. These histograms are strongly biased 
since they overweight the properties the most luminous 
galaxies in the samples. In order to correct for the 
Malmquist bias, so that histograms are weighted accord- 
ing to the true number density of galaxies, we have used 
the ratio between the LF derived in §[3] $(M), and the 
observed histogram of absolute magnitudes, 



AM 



,M,-M 
AM 



(15) 



As usual, the symbol 11 stands for the rectangle function, 

^^^^ { elsewhere. ^^^^ 

The index i in equation (|15p includes all observed galax- 
ies, whereas M and AM determine the centers and the 
widths of the histogram bins ^. We have to assume that 
each observed galaxy is actually a proxy for bi galax- 
ies, &i — 1 of which are not included in our observation 

''' Note that the histogram defined in equation (1151 difli'ers from 
that used in Fig. |4] because of the factor AM~^. This trivial re- 
scaUng is used here for convenience, allowing the corrected his- 
tograms to be normalized to the number density of galaxies. 



because they are too faint to exceed the apparent mag- 
nitude threshold. The magnitude of the galaxy causes 
the bias and, therefore, we assume that the bias function 
only depends on the absolute magnitude of the galaxy. 



6(M,). 



(17) 



This bias is precisely the reason why h{M) is not the LF, 
therefore, for AM small enough. 



$(M) 



AM ^ 



AM 



(18) 



Using equation ([TT]) . and assuming that 6(M) is a slowly 
varying function of M, 



6(M,)n( 



M, -M^ 
AM ' 



, ^^,Mi-M. , , 

Wn(^^), (19) 



and so one finds from equation (jlSp the expression for 
the bias function. 



b, = 6(M,) = ^{M,)lh{M,). 



(20) 



Let us denote as fc(p) the histogram of any parameter p 
constructed using the observed galaxies. 



(21) 



with A-p the bin-size and pi the value corresponding to 
the i-th galaxy. Then the bias-corrected histogram would 
be, 

^*(p)-^E^^n(^^""^ 



^ Ap ' 



(22) 



This definition guarantees that the corrected histogram 
of observed absolute magnitudes is the luminosity func- 
tion (cf. equations [TH] and [22] with p — M) and, there- 
fore, it is not difficult to show that all corrected his- 
tograms have the same normalization as the LF, i.e.. 



fc*(p) dp = Tio, 



(23) 



no being the number density of galaxies (see App. [X)l . 

We have applied equations ([20]) and ([22]) to restore the 
histograms in Fig. [4] The result is shown in Fig. [51 The 
means, standard deviations and modes of the new his- 
tograms are also included in Table [21 As expected, the 
mean absolute magnitudes of the restored histograms are 
much fainter than the observed ones. However, in addi- 
tion to this effect, the corrected histograms hint at in- 
trinsic metallicities significantly lower than the observed 
ones both for BCD and QBCD galaxies. The corrected 
QBCD galaxy colors are bluer than the colors of the ob- 
served set. These changes can be pin down to the correla- 
tions between luminosity, color, and metallicity existing 
in the original data set. Figure [9] shows how the QBCD 
galaxies tend to be less metallic as they become fainter, 
and a similar, but less marked trend, is also present in 
BCD galaxies. Such a relationship between metallic- 
ity and mass in dwarf galaxies (and so, between metal- 
licity and luminosity) was given by iPagel fc Edmuiidi 
( 1_981 ) , and later on by many others (see lKunth fc OstlinI 
2000i). Figure [9] includes the linear relationship found by 
Skillman et alT(il989, ) for nearby dwarf irregular galax- 
ies. Curiously enough, the slope is almost the same as we 
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(/ig> [mag arcsec-2] 12 + log(0/H) H„ EW [A] 

Fig. 8. — Histograms of magnitudes and colors for QBCDs (the solid lines) and BCDs (the dotted lines) once the Malmquist bias 
has been compensated for. The unrestored versions of these histograms are shown in Fig. |4| The labels in the ordinate axes describe 
the parameter that is represented. The histograms are normalized to the number density of galaxies. The dashed line in the plot of 
metallicities (12 + log[0/H] ) repro duces Fig. 3 in IKunth &: OstlinI 1120001) . which corresponds to the abundances of 425 H II galaxies 
cataloged bv lTerlevich et all 1119911) . Refer to Fig. |3] for further details. 
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Fig. 9. — Metallicity vs absolute Mg magnitude. The symbols 
stand for the mean values considering all the galaxies with similar 
Mg , whereas the error bars represent the standard deviation of the 
mean values. Both QBCD galaxies (the square symbols) and BCD 
galaxies (the asterisks) follow a trend so that the fainter the galaxy 
the more metal poor. The solid line is shown for reference and it 
corresponds to the law for nearby dwarf irregular galaxies found 
by Skill man et al. (1989 ). 



Fig. 10. — {/ig) — (fir) color vs absolute Mg magnitude. The 
symbols stand for the mean values considering all the galaxies with 
similar Alg , whereas the error bars represent the standard deviation 
to be expected for these mean values. The QBCD galaxies are bluer 
as they become fainter (the square symbols). The same trend is 
not so obvious for the BCD galaxies (the asterisks). 



find for QBCD galaxies. There is also a relationship be- 
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Fig. 11. — Metallicity vs color. There is a clear trend, which 
becomes steeper and sharper in the case of the QBCD galaxies 
(the square symbols). The solar metallicity has been marked for 
reference. Each symbol represents the average metallicity (12 + 
log[0/H]) considering all the galaxies with the same color. The 
error bars correspond to the standard deviation to be expected for 
these averages. 



tween color and luminosity, so that fainter QBCD galax- 
ies seem to be bluer; see Fig. [TU] Figures M and 
1101 actually reveal an underlying relationship between 
color and metallicity. It has been brough up in Fig. [Til 
showing how the redder the galaxy the larger the oxygen 
abundance. The relationship is more clear in the case of 
the QBCD galaxies (the square symbols), but it is also 
present among the redder BCD galaxies (the asterisks for 
(fig) — {^r) > —0.4). The relationships for QBCDs and 
BCDs meet at {fig) — (/i^) — —0.1, although they have 
different slopes. 

The number of known BCD galaxies with metallicities 
less than 1/20 the solar value is rather small (a dozen 
or so, according to iKunth fc Ostlin|[2000f) . Surprisingly 
enough, the list of 1609 BCD candidates we work out 
contains only three candidates below this metallicity. We 
interpret this result as a support of the metallicity esti- 
mate carried out in the paper, but it also reinforces the 
existence of a minimum H li region based metallicity (see 
IKunth fc Ostlinll200 0'. and references therein). 

The range of selected QBCD colors is rather broad and, 
therefore, our list of QBCD candidates seems to include 
all the range from early-type galaxies to late-type galax- 
ies. This fact can be appreciated in Fig. [T^ . which con- 
tains a color-color scatter plot. Our galaxies follow the 
sequence used in galaxy classification, with the reddest 
extreme corresponding to elliptical galaxies (early types), 
and the bluest extreme to irregular gala xies (late types) 
- com pare Fig. [T2k with, e.g.. Fig. 2 in iBershadv et ahl 
(12001 . One can distinguish two clusters or concentra- 
tions in Fig. [T^ (the dashed line arbitrarily separates 
the two sub-sets) . The existence of these two clusters or 
classes was to be expected, since they correspond to the 
bimodal distribution of colors found in t he large sam- 
ples of galaxies (e.g., iBalogh et al.l l200l : iBlanton et al.l 
|2005[ ). Hints of these two modes are also found as two 
small peaks in the histogram of surface brightness color 
shown in Fig. 2] (top right, the solid line, with the two 
peaks at (/ig) — {^r) — 0.4 and 0.7). Note how the BCD 
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Fig. 12. — (a) n—g vs g—r scatter plot for the QBCD candidates. 
The colors in the sequence span from those of early-type galaxies, 
to those of late-type galaxies. There seem to be two clusters in the 
distribution of points, which we have artificially separated by the 
dashed line, (b) Same as (a) but showing the BCD candidates. 

candidates occupy the bluest extreme of the color-color 
plot; cf. Figs. flA and fT2b. 

Most of the QBCD candidates show Ha in emission, 
meaning that even if the star-formation is reduced with 
respect to the BCD candidates, it is not absent. Only 
some 3200 QBCD galaxies out of the 21493 candidates 
show Ha in absorption. They have (negative) equiva- 
lent widths of a few A and, therefore, these objects are 
included in the bin of zero equivalent width in Fig. |4l 
bottom right. Using the Ha flux as a proxy for star for- 
mation (jKennicu ttl ll998i) . our QBCD candidates turn out 
to be one order of magnitude less active in forming stars 
than the BCDs. This estimate is worked out in § 15. H 
equations ^ and 

5. DISCUSION 

So far we have described our search avoiding interpret- 
ing the results. This section, however, is fully devoted to 
interpretation and so, admittedly speculative. We will 
discuss how the results in previous sections are consis- 
tent with the BCD candidates and the QBCD candidates 
being two different phases in the same sequence of galac- 
tic evolution. Using as a working hypothesis that BCDs 
change into QBCDs, and vice versa, we examine obvi- 
ous flaws and constraints. Showing the consistency of an 
hypothesis does not prove it to be correct - other alter- 
natives cannot be discarded. It just shows the working 
hypothesis to be viable. Three main results point out in 
this direction. First, we choose the QBCDs so that their 
colors, concentration indexes and luminosities are equal 
to those of the galaxies underlying the BCDs. Second, 
the LF of QBCDs is conformable to the LF of BCDs 
(§ O. Third, the BCD galaxies turn out to be QBCD 
galaxies with the largest (specific) star formation rates 

(§E2l). 

Matching LFs forces the QBCD phase to last some 
30 times longer than the BCD phase (§ 13. ip . Since the 
BCD starbursts are so short-lived, there should be many 
episodes of BCD phase during a galaxy lifetime. As- 
suming a single 10 My long starburst per BCD phase, 
there should be one of such BCD phases each 0.3 Gy. 
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Longer BCD phases are in principle possible by concate- 
nation of starbursts, but population synthesis modeling 
of BCD spectra does not favor ex tended periods of star- 
formation fe.g.. lMas-Hesse fc Ku nth 1999). If BCDs are 
transformed into QBCDs and vice versa, each QBCD has 
experienced several BCD phases because the other alter- 
native does not seem to be consistent with the results of 
the search. Let us examine this alternative possibility, 
namely, that the BCD phase represents the only episode 
of intense star formation during the galaxy lifetime. Sta- 
tistically, all QBCD galaxies that we detect today have 
sufferred a BCD episode during the last 270 My (equa- 
tions [n] and [M], with tbcd =10 My). If this is the 
only such episode per galaxy, there should be no BCDs 
among galaxies observed at lookback times larger than 
270 My, or at redshifts > 0.02 (= 270 My/Ho). How- 
ever, a large fraction of the sample of BCDs selected in 
ij 12.21 has redshifts larger than this limit. We are forced 
to conclude that the same host galaxy undergoes several 
BCD phases. 

Although the recursive BCD phase scenario explains 
a number of important observables, other properties of 
QBCDs and BCDs do not seem fit in the picture so well. 
The purpose of the rest of the section is to discuss (and 
hopefully clarify) some of the most obvious difficulties 
posed by the scenario. 

5.1. Hi consumption timescale 

One may think that the repeated starbursting of a 
dwarf galaxy quickly exhausts the gas reservoir, thereby 
making the whole picture inconsistent. However, the low 
surface brightness galaxies like the QBCDs are gas-rich 
(jStavelev-Smith et al.l[l99l iKunth fc Osnh] [20(X^ with 
fuel to power the star formation during several Hubble 
times. This claim can be supported by working out the 
H I consumption timescale, thi, which is defined as the 
time required to transform into stars the neutral hydro- 
gen of the galaxy. A/hi, i-e.. 



THI = Mhi/SFR, 



(24) 



with SFR the star formation rate. On the one hand, 
dwarf low surface brightness galaxi es have about one so- 
lar mass of H i per solar luminosity (jStavelev-Smith et al.l 
fl992i) . exphcitly. 



10g(AfHl/Mo) = -0.4 {Mg - MgQ), 



(25) 



Mq being one solar mass, and A/g© the solar absolute 
luminosity in the g filter. On the other h and, the SFR 
scales with the Ha luminosity (e.g. jKenni cutt 1998), i.e., 
with the product of the H^ equivalent width times the 
galaxy luminosity. Since the observed BCD Ha equiv- 
alent width is uncorrelated with the galaxy luminosity, 
the BCD SFR scales with the luminosity, that is to say, 

log(SFRBCD/SFRi3CD) = -0.4(Afg - AQ, (26) 

where SFRg^j-, stands for the SFR of a BCD galaxy of 
magnitude M'g. In our two-phase scenario, the galaxy has 
short periods of star formation lasting tbcd , interleaved 
with long periods of quiescence lasting tqbcd (§ 13. ip . 
Consequently, the effective time-average SFR is. 



g-pj^ ^ Tbcd SFRbcd + '^qbcd SFRqbcd ^^^^ 

TBCD + TQBCD 



SFRBCD(^^+e), 



Tqbcd 



with 



SFR, 



QBCD 



SFRbcd 



(28) 



As judged from the ratio of Ha equivalent widths in Ta- 
ble[2l and the difference of luminosity between BCDs and 
QBCDs (equation P), the SFR during the BCD phase 
is some ten times larger than during the QBCD phase. 



I.e., 



10. 



(29) 



By combining equations and ^7^, one 

finds the consumption timescale to be independent of 
the galaxy luminosity. 



logTHi ^ -O.4(A/'-Afg0)-log 



SFRl 



BCD 



Ma 



' Tbcd 
'Tqbcd 



(30) 

Although BCD star formation bursts are intense for a 
dwarf galaxy, the SFRs of BCDs are rather modest, i.e., 
less than lA ^^y"^, and typically one order of mag nitude 
smaller fe.g. lSage et al.lll992l: [Mas-Hesse fc Kunthlll999t) . 
Taking SFRbcd < l-^oY"^ for the brightest galaxies of 
our sample, M'^ ~ —19, and using the ratio of timescales 
between the phases as provided by equations (|lip and 
(mD, equations (USD and ^ yield. 



THI > 30 Gy. 



(31) 



Consequently, the hydrogen existing in QBCD galaxies 
allows the sequence BCD-QBCD to last for a few Hubble 
times (= I /Ha =14 Gy). 

A final comment is in order. The scenario of a recur- 
sive starbursting is not in conflict with the commonly 
accepted view of a SFR decreasing du ring the last 10 Gy 
(|Madau et al.lll996t iLiTlv et al.lll996( l. This drop refers 
to massive galaxies. However, the kind of dwarf galax- 
ies in our QBCD sample, with stellar masses less than 
10^° A/q, has maintained a SFR either stationary or 
slightly increasing with time (see iHeavens" "eFaI1[200l 
Fig. 1). 

5.2. Metallicities 

The metallicities of our QBCD galaxies are system- 
atically larger than the metallicities of the BCD galax- 
ies. This result seems to be in conflict with the recursive 
BCD phase scenario, since the BCD starburst starts off 
with a metallicity lower than the metallicity of its host 
galaxy (see Tabl£[2|). Although this disagreement may 
indicate a real flaw in the overall picture, one can also 
think of various ways to circumvent the apparent incon- 
sistency. For example, the BCD episodes may involve 
fresh low metallicity gas accreted by the host galaxy 
during the periods of quiescence. Note that the BCD 
episodes are very conspicuous from the point of view of 
the luminosity, but they involve moderate masses, which 
hardly exceed 10^ Mq (e.g. Mas- Hesse fc Kunthi ri999). 
A sustained infall such as those postulated in the lit- 
erature to solve various problems in galaxy evolution 
IMq pc-2 Gy"^ see iDalcantonl l2007l and refer- 
ences therein) provides 10^ Mq in only 50 My, even for 
a small 2.5 kpc galaxy. Masses of 10^ Mq are also typi- 
cal of the high velocity clouds hounding our galaxy (e.g.. 
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iWakker et al1[200l . If these gas reservoirs are common 
in QBCDs, they may develop major star formation events 
when merging with the galaxy. 

The fresh gas glowing during the BCD phase does 
not necessarily have to come from outside. Metal-poor 
gas may exist in place in the galaxy (§ 15. ip . Then 
the oxygen metallicity assigned here to QBCD galax- 
ies may not reflect the metallicity of this gas, but the 
metallicity of a few metal-polluted H ii regions remain- 
ing from the last starburst. It is relatively easy to en- 
rich with metals the small fraction of galactic gas un- 
dergoing starbursts. For example, the metallicity of a 
pure gas cloud increases from the level observed in BCDs 
(12 -I- \og{0/H) ~ 8.24; Table [2]) to the level in QBCDs 
(~ 8.61), if 60% of the original mass is transformed into 
starts.* The metal-polluted gas produces emission lines 
during the extended recombinatio n phase of the H ii re- 
gion (e.g.. iBeltrametti et al.|[l98^ . and it may also give 
rise to a secondary generation of stars. The light from 
these aging H ii regions renders high metallicity mea- 
surements. Eventually, the star-processed gas mingles 
with the metal-poor gas of the galaxy, but this mixing 
leaves the original metallicity almost unchanged^. We 
note that this explanation may be in conflict with the IR 
excess detected in the halos of Blue Compact Galaxies 
by B ergvall fc OstlinI (|2002t l. These halos can be identi- 
fied with the BCD host galaxies and, therefore, with the 
QBCDs. The colors cannot be explained with a normal 
metal-poor stellar population like the Milky Way halo, 
but an excess of lo w-mass stars with significant metallic- 
ity is required (see IZackrisson et al.l 12006) . The nature 
of the disagreement requires further investigation, and it 
may be due to the fact that these Blue Compact Galaxies 
with IR excesses are not dwarfs according to the criteria 
used in our work. 

There is yet another possibility to reconcile the mctal- 
licities of the two galaxy types. Only those QBCD candi- 
dates with metallicity similar to that of the BCD galax- 
ies would be part of the BCD - QBCD sequence. The 
QBCDs with metallicity like the BCDs represent 20% 
of the sample (see Fig. U) and, therefore, the timescale 
of quiescence must be shortened by a factor of five to 
comply with equation (|14p . The same happens with the 
H I consumption timescale worked out in equation (|3ip . 
These QBCD galaxies may keep their low metallicity lev- 
els because the metals are expelled by ga lactic winds 
(iHeckmanI 120021: If enorio-Taalc ct a ll l2003t see, how - 
ever. ISihch fc Tenonio-Tagle,2001; .Legrand et al.ll200lD . 
because they are recycled into new stellar generations 
within superstar clusters, w ithout significant mixing with 
the galactic gas (jTenorioJCagle et al. 2005), or simply be- 
cause of the dilution with pristine gas described above. 

5.3. Starburst triggering 

We select isolated QBCD galaxies, therefore, the star- 
burst triggering mechanism to turn them into BCDs 
cannot be galaxy-galaxy interactions, mergers or harass- 
ment. How, then, are the BCD bursts triggered? Sev- 

® This estimate assumes the closed box evolution of a purely 
gaseous cloud (e.g. Tinslev 1980) with the standard oxygen yield 
(~0.003; see, e.g., Pilvu gin et al.n2004i ). 

® If 1% of gas with QBCD metallicity is mixed up with 99% of 
gas with BCD metallicity, then the resulting change of metallicity 
is only A[12 -|- \og{0/H)\ ~ 0.006, i.e., insignificant. 



eral possibilities are available. If the fresh gas falls on the 
QBCD galaxy over extended periods of time, intense star 
formation cannot be triggered until the gas density ex- 
ceeds the requ ired threshold (the so-called Schmidt law; 
ISchmidtl 119 59*). Waiting for enough gas to accumulate 
could explain periods of latency. If the gas falls in during 
short episodes (like the collision with a large gas cloud), 
the hitting of the gas itself may induce star formation. 
Then the QBCD timescale would be given by the char- 
acteristic time between cloud collisions. If rather than 
coming from outside, the fresh gas is part of the galaxy, 
internal galactic structures like bars may periodically ex- 
cite star formation. In the case of a gas-rich galaxy the 
triggering may also be due to perturbations of the galac- 
tic material by dark matter clumps. These elusive struc- 
tures are predicted in lots by the cold dark matter sim - 
ulations of galaxy formation (e.g. jDiemand et al.ll2007l ). 
and they are expec ted to disturb the quiet evolution of 
galactic disks (e.g.. iKazantzidis et al.l[2007[ ). 

5.4. QBCD characterization criteria 

The criteria used to characteriz e QBCDs com e from 
the B CD host galaxies analyzed bv lAmorm et al.l (|2007l . 
[2008I) . but we may have used different criteria taken from 
other works. The question arises as to whether the selec- 
tion of candidates critically depends on this assumption. 

The BCD host galaxy properties are used to con- 
strain colors and Sersic indexes (see ij 12.11 and Table [2]). 
The lower limit luminosity, also taken from the BCD 
host galaxies, turns out to be unimportant since the 
small number low luminosity objects prevents any se- 
rious influence of this criterion on the selection (see 
Fig. [4]). The range of visible g — r colors that we use 
agrees with the colors that dif ferent studies have as- 
signed to the hosts of BCDs (e.g., IPapaderos et al.lll996bl : 
iGil de Paz fc Madord [20051 ) ■ even in those cases where 
a significant infra-red excess has been detected (e.g., 
iBergvall fc Ostlh]|[200l . In this sense the color crite- 
rion does not seem to be questionable and, therefore, 
any work would have rendered similar QBCD candidates. 
The selecti on of Ser sic indexes is niore influential. Stud- 
ies like BcrgvalL&Ostlin (200l) or lGil de Paz fc Madord 
|2Q05)) find indexes as large as n =10-20, whereas we take 
n < 3 (Table [2]). However, there are two good reasons to 
prefer low— n values. First and most important, the ob- 
ser yed dwarf ga laxies present low— n values (e.g. . Fig. 11 
in ICaon et al.l 120051 : iGraham fc Guzmal] l2003f ). and if 
we search for galaxies that may be Blue- Compact- _Du;ar/ 
galaxies during quiescence, they must be dwarfs too (see 
equation [T].) Second, the large n values are obtained 
from ID fits to azimuthally averaged luminosity profiles, 
and the n thus obtained depends critically on the range 
of radii used to characterize the galactic outskirts (e.g., 
iCairos et al.ll2003l: IGil de Paz fc Madordr2005 ). The 2D 
fits providing low indexes ar e fairly more robus t and, 
therefore, to be preferred fsee lAmorm et al.ll2007l ). 

6. CONCLUSIONS 

The starburst characteristic of BCD galaxies does not 
last for long. Once this episode is over, the BCD galaxies 
should dim their surface brightness to become quiescent 
BCD galaxies (or QBCD galaxies). Although QBCD 
galaxies are to be expected, they have not been iden- 
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tified yet (see § [Ij . The present work describes an ef- 
fort to find them among the galaxies in the SDSS/DR6 
database. The properties of the QBCDs have been taken 
from the sample of BCD host galaxies characterized by 
lAmorm et alflpOOTl . [200 8) (§EI|). We find 21493 QBCD 
candidates, therefore, QBCDs seem to be fairly common 
in the local universe. In order to have a proper refer- 
ence to compare with, a complete sample of BCD galax- 
ies was selected too. It comprises 1609 BCD candidates 
(§ 12. 2p . Since the two samples were selected from the 
same database using analogous criteria, the comparison 
between the two sets is relatively free from the bias that 
our selection may have. 

We compute Luminosity Functions (LFs) for the two 
samples (§[3]). In addition, we estimate their main prop- 
erties before and after correcting for Malmquist bias. 
These properties can be summarized as follows, 

• There are around 30 QBCD candidates per BCD 
candidate. We infer this ratio by comparison of the 
LFs for QBCDs and BCDs. The two LFs are very 
similar, except for the global scaling factor, and the 
expected 0.5 magnitude dimming 13. ip . 

• The surface brightness of the QBCDs is typically 
one magnitude fainter than the surface brightness 
of BCDs (Table E]). 

• QBCD candidates are, on average, 0.4 magnitudes 
redder than the BCD (Table H {^J,g) - (/z^)). 

• QBCD candidates have an H ii region based oxygen 
metallicity 0.4 dex higher than the BCD candidates 
(Tabled'. 

• The QBCD metallicity increases with the luminos- 
ity, following the well known trend for dwarf galax- 
ies. 

• The QBCD metallicity also increases with the 
color, so that the redder the galaxy the larger the 
measured metallicity. 

• 75 % of the BCD candidates are also part of the 
QBCD sample (§1221). Roughly speaking, the BCD 
sample represents the fraction QBCD galaxies hav- 
ing the largest specific SFR (SFR per unit of lumi- 
nosity) . 

• There are around three dwarf galaxies per QBCD 
candidate (§[3|), which renders one BCD galaxy ev- 
ery ninety dwarf galaxies. 

The overlap between BCD galaxies and QBCD galax- 
ies is consistent with the two sets forming a single con- 
tinuous sequence, with the most active QBCD galaxies 
being BCD galaxies. The agreemet between their LF 
shapes, and the ratio of number their densities, support 
the commonly accepted view that BCD galaxies undergo 
short bursts of star formation separated by long quies- 
cent epochs. However, the fact that the QBCD metallic- 
ity is higher than the BCD metallicity poses a problem to 
such an episodic starbursting scenario, which we try to 
circumvent with various plausible explanations in § 15.21 

The repeated starbursting scenario predicts a number 
of independent observables whose testing is important 



but goes beyond the scope of the paper. They represent 
a natural extension of the present work. In order to il- 
lustrate the possibilities, we will outline two examples. 
QBCDs must have a stellar population corresponding to 
short star formation episodes in between quiescent gaps. 
In principle, one can distinguish between a (low level of) 
continuous star formation and a mor e violent episodic 
star formation with a period of 0.3 Gy. iLeonardi fc Ros j 
(jT996) put forward a spectroscopic index to detect post- 
starburst galaxies. The technique is well suited for de- 
termining the time elapsed from t he last (young) star- 
burst, and it has been calibrated bv ILeonardi fc WorthevI 
(|2000l ). Depending on the noise level, one can apply the 
method to selected QBCD SDSS spectra, or to averages 
of similar QBCD spectra. Another testable prediction of 
the scenario has to do with the metallicity inferred from 
emission lines in QBCD galaxies 15. 2p . It should over- 
estimates the true galaxy metallicity and, in particular, 
the metallicity of the stellar content. Studies of stellar 
metallicity can be carried out using integrated galaxy 
spectra, provided that they have enough signal-to-noise 
ratio and spectral resolution to show absorpt ion lines 
(e.g. iTerlevich et al.lfl990t IWorthev et al.lll994f ). These 
two testable predictions provide a flavor for other tests 
to come. 
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APPENDIX 

LUMINOSITY FUNCTION ESTIMATE 

The LP, $(M), is defined as the number of galaxies with absolute magnitude M per unit of magnitude and unit of 
volume. We use it to describe the local universe and, therefore, <I>(M) is not expected to change with the position in 
space. The number of galaxies per unit volume, no, is just, 

710 = / $(M)dAf, (Al) 
Jmi 

where only absolute magnitudes in between the limits A// < A4 < Mu are considered. Using the definition of LP, one 
can easily write down the number of galaxies to be expected in an apparent magnitude limitted catalog, namely, 

/ $(M)Kiax(Af)dM. (A2) 

J Ml 

As usual, the symbol Vmax(-M) stands for the volume of universe covered by the catalog where galax ies of absolute 
magnitude M have apparent magnitudes smaller than the catalog threshold fe.g. . iTakeuchi et al.l 200dl). 



Maximum likelyhood estimate s of LPs are fav ored in the current literature (e.g.. lLin et al.lll996HBlanton 
They were pioneered by Efstath iou et al.l (|1988f l and, among the quoted advantages, they present a number of desirable 
asymptotic error properties'^. Note, however, that these methods are not unbiased (the expected value of the estimate 
is not necessarily the parameters to be estimated). The method is based on the following principle; given the redshift 
C of an observed galaxy, the probability that it has an absolute magnitude M is, 

, /-M„„(C) 

P(M I C) = «'(M) / / $(M') dM', (A3) 

where the normalization gives all the galaxies that we are allowed to observe at C. Since our catalog is limited in 
apparent magnitude mmax, 

m < TOinax, (A4) 

only galaxies bright enough would be observable at this redshift. Considering the definition of distance modulus, 

DM{C) = m- M, (A5) 

then, 

M < M„,ax(C) — 772 max DM{C). (A6) 

Similarly, the sample would have a minimum apparent magnitude TOmin (e.g., given by the brightest galaxy in the 
sample), which sets the minimum magnitude to be observed, 

M>Af^i„(C)= TOmin -CM(C). (A7) 

Assuming that the probability of observing each galaxy is independent from the rest of galaxies, the likelyhood function 
£ is just the product of the probability of observing each galaxy, 

C = njP{M,\Q, (A8) 

where the index j spans from 1 to N. The maximum likelyhood estimate maximizes C, which is equivalent to 
maximizing its logarithm. 



log / <^>{M')dM' . (A9) 

The next step consist in parameterizing the LP in terms of a number K of free parameters Xi, 

<^{M)^4>{M,Xi), (AlO) 

with i = \,...K. V arious representations 0(Af , Xj) can be found in the literature, e.g., a step wise function 
(jEfstathiou et al.| 119881 ). a Schechter function ()Lin et al.lll996( ). a collection of Gaussians (jBlanton et al. 2003b). We 
choose yet another representation, namely, natural cubic splines. The reasons are (1) the LP is simple and fast to 
compute, (2) the integral of the LP is also fast to compute since it follow directly from the spline interpolation, and 
(3) it automatically provides a smooth '^{M). Speed is always an appealing feature since the maximization of log£ is 
carried out iteratively, and requires many evaluations of the LP. 

The normalization of the LP is not constrained by the likelyhood function, which is independent of a global scaling 
factor (equation |A9j ). Another complementary method is required to est i mate the number density of galaxies uq. 
Such method is often the minimum variance estimate by iDavis fc Huchral ()1982[ ) . We use a simple version of such 

They are consistent, i.e., they tend to the parame ter to be est imated when the sample increases, and the distribution becomes a 
normal of minimum variance for large samples (see, e.g., IMartinlllQTll . § 7.2 ). 
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estimate, where all the galaxies in the sample are equally weighted, thus avoi ding assuming a particu lar covariance of 
the galaxy catalog. I t corresponds to the so -called ns in the original paper by lDavis &: HuchrU ()1982D . and it has been 
used elsewhere fe.g.. iBolzonella et al.|[2002l ). The expression can be derived by combining equations (jAl[) and (|A2[) . 
which yield, 

uo^n/ [ \{M)V^UM)dM. (All) 

' J Ml 

All items in the right-hand-side of the previous expression are known; the normalized LF ip{M), 

ip{M)^^{M)l I <^{M')dM\ (A12) 
' Jmi 

is provided by the maxim um likelyhood proced ure, whereas Vmax(Af) follows from the magnitude limit and the solid 

angle of the catalog (e.g-. TTakeuchi et al.ll2000h . 

The actual max imization of the likelyhood in equation (|A9[) is carried out using the standard Powell method 
(jPress et al.lll988[ ). Error bars are assigned by bootstrapping (e.g.. iMoore et al.]|2003t [Blanton et al.ll2003bl ). where 
one constructs bootstrap re-samples by choosing at random galaxies from the original data set. Then the application 
of the LF retrieval procedure to all bootstrap re-sa mples yields a set of LFs with the spread of values to be expected 
from the true error distribution (jMoore "eralll2003l) . We use the standard deviation of such bootstrap distribution as 
our error bars. 

Numerical tests have been carried out to check the procedure. We choose a large number of galaxies (10^-10^) with 
random absolute magnitude according to Gaussian or uniform distributions. These galaxies are randomly uniformly 
spread in space within a sphere of radius C = 0.35. Apparent magnitudes are computed using these redshifts, and 
galaxies fainter than the assumed catalog cutoff are drooped from the sample (we take rrimax = 17.7). This biased 
sample is then used to fed the procedure. Two examples of true and restored LFs are given in Fig. [T31 The solid 
lines correspond to our maximum likelyhood estimate, whereas the dashed lines show the true histogram derived from 
the synthetic data before the apparent magnitude threshold is introduced. The two of them agree within error bars, 
which are given as shaded areas. M oreover, the ma.ximun likelyhood estimates also agree with the 1/Vmax estimates 
included in the same plot (see, e.g.. iTakeuchi et al]|2000l for a description of traditional 1/Vmax method bv lSchmi"dtl 
Il968| ). The agreement is not specific of these particular realizations but is a general property, and it seems to hold 
both independently of the original distribution, and the number of galaxies in the sample. In order to illustrate the 
magnitude of the correction carried out by the LF retrieval routine, Fig.[T3]also includes the distributions of magnitudes 
of the some 10"* galaxies given to the program (the dotted lines). 
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Fig. 13. — Numerical tests to our maximum likelyhood LF estimate. The various types of lines and symbols correspond to the true 
LF (the dashed line), the restored 1/Vmax LF (square symbols), and the maximum likelyhood LF (the solid line, with the shaded region 
representing the error bars). The dotted line corresponds to the original data offered to the LF retrieval routine, (a) Case where the true 
LF is a uniform distribution between magnitudes -23 and -17. (b) Same as (a) except that the underlying LF is a truncated Gaussian. 
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